The hierarchy of multiple many-body interaction scales in high-temperature 

superconductors 
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To date, angle-resolved photoemission spectroscopy has been successful in identifying energy 
scales of the many-body interactions in correlated materials, focused on binding energies of up 
to a few hundred meV below the Fermi energy. Here, at higher energy scale, we present improved 
experimental data from four families of high- Tc superconductors over a wide doping range that 
reveal a hierarchy of many-body interaction scales focused on: the low energy anomaly ("kink") of 
0.03-0.09eV, a high energy anomaly of 0.3-0.5eV, and an anomalous enhancement of the width of 
the LDA-based Cu02 band extending to energies of « 2 eV. Besides their universal behavior over 
the families, we find that all of these three dispersion anomalies also show clear doping dependence 
over the doping range presented. 
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I. INTRODUCTION 

Many-body interaction is a key to understanding novel 
properties of quantum matter. As an extreme example, 
the complexity due to charge, spin, and lattice inter- 
actions in high- Tc superconductors makes it difficult to 
identify the essential microscopic ingredients for the basic 
model - a reason behind the current debate on the mech- 
anism. The energy-momentum dispersion relationship 
measured by angle-resolved photoemission spectroscopy 
(ARPES) provides an excellent tool for identifying these 
scales. Energy scales where these interactions are man- 
ifest usually provide important insights into the nature 
of the interactions. At an early stage, APRES proved 
to be successful in identifying the energy scale of the 
d-wave gap in high-Tc superconductorsii^. To date, a 
focus of the discussion in ARPES has been on the na- 
ture of electron-boson coupling, which manifests itself in 
the form of the low energy anomaly ("kink") near 0.03- 
O.OQeV^iii^i^ii^i^iia. However, little attention has been 
paid to the features at higher binding energies. 

Given the strong many-body interactions and com- 
plex band structure, one expects only incoherent fea- 
tures and complex spectral weight modulation at high 
energy. The fact that one can see neatly defined mo- 
mentum dependent features which are robust against 
doping and measuring conditions (e.g. photon en- 
ergy) is unexpected and thus provides a new oppor- 



tunity to understand many-body effects beyond what 
have traditionally been the points of focus by ARPES, 
other than some early work on oxygen p bands^*'*^. 
Here, we report ARPES experiments on four fami- 
lies of high-Tc cuprates over a wide range of dop- 
ings: Bi2Sr2Cu06 (Bi2201), Bi2Sr2CaCu208 (Bi2212), 
Ba2Ca3Cu408(05,Fi_5)2 (F0234), and La2-r.Sr,Cu04 
(LSCO). All of their measured energy-momentum dis- 
persion relationships reveal the simultaneous presence of 
three energy scales, marked as 1-3 in Fig. Id: 1) the 
band bottom, at nearly 2 eV in optimally-doped Bi2201, 
which is deeper than predicted by LDA calculation, 2) 
the high-energy anomaly (HEA) at ~ 0.3-0.5 eV (green 
arrow) and 3) the low-energy "kink" (LEK) near 70 meV 
(red arrow) which can be better seen in an enlarged en- 
ergy window. Various aspects of these energy scales, in- 
cluding the peculiar doping dependent effects, will be dis- 
cussed to obtain insights on many-body interactions and 
their interplay in cuprates. 



II. EXPERIMENT AND LDA CALCULATION 

We have measured four families of high-Tc cuprates: 
Bi2Sr2Cu06 (Bi2201), Bi2Sr2CaCu208 (Bi2212), 
Ba2Ca3Cu408(05,Fi_5)2 (F0234>M, and La2-:rSr:rCu04 
(LSCO). The Bi2201 samples are optimally-doped of 
Tc = 35K and non-superconducting overdoped. The 
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FIG. 1: (a) Raw ARPES spectrum, along (0, 0) to (tt, tt) direction, of non-superconducting overdoped Bi2201 at T = 30 K 
while (b) and (c) represent its MDCs and EDCs respectively, (d) the raw data is normalized with angle-integrated EDC profile 
for clearer view and the numbers 1-3 mark the wider band width, the high and low-energy anomalies respectively. Normalized 
ARPES spectra are compared with LDA calculations as follows: (e) OP Bi2201 with = 35K (T=45K), (f) OD Bi2201, 
non-superconducting (T=30K), (g) OP Bi2212 with = 96K (T = llOK, LDA from Ref^) and (h) OD Bi2212 with Tc = 
65K (T=76K). Inset in Fig. la shows the momentum space of the data. Note that the LDA bands in panel (h) are obtained 
by rigidly shifting the bands in (g) to account for the correct doping level. It is not clear if feature B should be matched to the 
top of the band structure in (e) or (f) at F; however, this uncertainty is not important for our argument which only requires a 
relative shift in going from the OP to the OD case. 



Bi2212 samples are optimally-doped of To = 92K and 
overdoped of To = 65K. The F0234 samples are of To 
— 60K. And, the LSCO samples has a wide range of 
dopings: X = 0.03, 0.05, 0.063, 0.07, 0.075, 0.09, 0.12, 
0.15, 0.22 and 0.3. The measurements were carried 
out on beamline 10.0.1 at the ALS, using a Scienta 
R4000 electron energy analyzer. This analyzer has the 
advantage of a large-angle window which can cover the 
band dispersion across the Brillouin zone as shown in 
Fig. 1. We stress that the wide angle scan allowed 
us to record above data without resorting to manual 
symmetrization. The photon energies are 37, 40, 41,.. .,45 
and 55 eV. The energy resolution between 12 and 20 
meV was used for various measurements on different 
samples, and the angular resolution is 0.3 degree. The 
samples were cleaved in situ in vacuum with a base 
pressure better than 4 x 10^^^ torr. The samples were 
measured both in normal and superconducting states. 



III. RESULTS 



Fig. la shows the raw ARPES image of the strongly 
overdoped Bi2201 sample while its raw momentum- 
distribution-curves (MDC) and energy-distribution- 
curves (EDC) are shown in Fig. lb and c, respectively. 
To see the band near the bottom more clearly, the raw 
ARPES image is divided by its profile of angle-integrated 
EDClS., as shown in Fig. Id; this procedure will not 
change the MDC-peak position at any given energy. We 
also note that this renormalization procedure which is 
used in Fig. Id-h and Fig. 2 is only for the purpose of 
displaying the lower and higher energy features together 
since otherwise the intensity at the higher-energy region 
will be too high to have a reasonable displaying contrast. 
The raw data without this renormalization procedure can 
be seen in Fig. la-c and Fig. 6a-i. 



LDA results here are based on full-potential well- 
converged computations for the appropriate lattice struc- 
tures, described in greater details in Ref.l^ andJ^. 



As marked by numbers 1-3 in Fig. Id, we will focus 
on how our new data reveal the simultaneous presence of 
three energy scales as follows. 
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FIG. 2: The comparison of normalized ARPES spectra and 
LDA calculations along (a) (tt, tt) to (0,0) or nodal direction 
and (b) (0,0) to (tt, 0) or antinodal direction in OD non- 
superconducting Bi2201 system, (c) shows the agreement of 
ARPES spectrum and LDA calculation of high-energy band 
when shifting the LDA by 450 meV to the higher energy, (d) 
The comparison of the spectra at photon energies, E^ = 40, 45 
and 55 eV where the green shaded region denotes the energy 
scale of HEA and blue shaded region indicates the top of the 
band B. We note that for E,y = 55 eV, due to the matrix 
element effect, the intensity of the left band B is large and 
therefore we have adjusted the color scale so that we can see 
the energy scales of both band A and B clearly. 



A. anomalous enhancement of the LDA-based 
Cu02 band width 

As shown in Fig. le-lh, ARPES spectra of optimally- 
doped (OP) and overdoped (OD) samples of Bi2201 and 
Bi2212 systems are overlaid on the corresponding LDA 
calculations. As seen in all measured samples, the first 
peculiar feature, especially at lower doping, is that the 
ARPES band width is found to be wider than LDA calcu- 
lation. This is anomalous as one expects interactions to 
enhance the mass and reduce the band width. By extrap- 
olating the band (e.g. the red dashed line in Fig. Id), one 
can get an estimate of the band width that is suitable for 
qualitative discussion. With doping, the discrepancy be- 
tween the band widths obtained from ARPES and LDA 
seems to be reduced. 

Additional evidence that these high energy dispersions 
still contain useful information comes from band B in Fig. 
la which shows a maximum at F point near 1 eV. This 
band has a correspondence to an LDA band and thus 



FIG. 3: (a) shows MDC-peak dispersion plotted on top of 
ARPES spectrum and (b) shows corresponding MDC width 
in OP Bi2201 system (Tc = 35K). (c)-(f) show the MDC- 
derived dispersions of Bi2201 (Tc = 35K), Bi2212 (Tc=65K), 
LSCO(Tc=38K) and F0234 (Tc=60K) respectively. Addi- 
tionally, the temperature dependence of Bi2212 and LSCO 
dispersion is shown in (d) and (f). 



provides confidence in the data at higher energy scales, 
which has been largely unexplored in the cuprates. From 
the LDA calculation with orthorhombic distortion, this 
band is the band at Y point (left arrow. Fig. 2c), which 
is folded around the (7r/2,7r/2) point. We then compare 
the LDA and ARPES top part of this band B at F point 
by shifting the LDA band down. A good agreement of 
ARPES and LDA of this concave-down band (see Fig. 
2c) can be obtained if the LDA is shifted down « 0.45 
eV for OD Bi2201 sample and the shifted energy increases 
to « 0.8 eV for OP Bi2201 samplc^^ leading to a filled 
band width near 2 eV. A similar behavior is also observed 
in the Bi2212 systen*ii (Fig. Ig-h). This band width en- 
hancement was also seen earlier in undoped Ca2Cu02Cl2 
(CCOC) as its high-energy dispersion matches with the 
LDA calculation shifted by 0.7 eViS. 



B. high energy anomaly of 0.3-0.5 eV 

Next, we discuss the high energy anomaly (HEA) near 
0.3-0.5 eV. We extract the MDC peak position by fit- 
ting to Lorentzian curves, as shown by the red and blue 
curves in Fig. 2a and 2b. We note that especially at 
high-binding energy, MDC peak position may not rep- 
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FIG. 4: The doping dependence of high-energy dispersion (red circles) in LSCO system at T = 20K at doping x = 0.03 to 0.30. 
LDA curves for various dopings (blue lines) are obtained by appropriate rigid shifts of the computations for x=0. 



resent the real dispersion. Explained in appendix, a full 
2D analysis, which directly extracts the spectral function, 
A{k,u}) and matrix element term separately at once, can 
avoid the problem of the MDC or EDC analysis alone. 
However, since the kink-like structure is so large, the 
MDC-derived dispersion should be able to approximately 
identify the energy scale of the HEA. To check the 3 di- 
mensional behavior of the band A and B, we perform the 
measurements at various photon energies 40, 41, 45 
eV and 55 eV, probing different perpendicular momenta 
kz- As shown in Figs. 2d, we find that HEA scale and the 
top of Band B are not very sensitive to the photo energies 
while as well these energy scales in our LDA calculations 
do not show strong kz dependence. 

This HEA is present in various cuprate families. While 
earlier seen in undopcd CCOC^*, the MDC-pcak-position 
graphs of Bi2201, Bi2212, LSCO and F0234 plotted 
in Fig. 3, reveal its universality. The energy scales 
are around 0.3-0.4 eV in Bi2201, Bi2212 and LSCO 
while around 0.5-0.6 eV in F0234 (for electron-doped 
band). HEA persists in both superconducting and non- 
superconducting samples, albeit its strength depends on 
doping. Fig. 4 shows the plots of MDC-peak position of 
LSCO samples which cover a wide doping range, x = 0.03 
to 0.30. From the figure, the HEA energy does not change 
much with doping. However, if we define the size of 
HEA to be the difference between the MDC-derived and 
LDA dispersion, it increases upon doping in this range. 
Similar doping-dependent behavior is also observed in 
Bi2201 and Bi2212 samples which cover a narrower range 
of doping. For superconducting samples, the HEA per- 
sists above and below Tc. consistent with ARPES data, 
from in-plane optical conductivity of Bi2212^^, Norman 
and Chubukov recently report that the real part of the 
self-energy is large with a maximum value around 0.3-0.4 



eV^. 



C. low energy "kink" of 0.03-0.09eV 

Finally, the low energy "kink" (LEK) around 0.03-0.09 
eV is indicated with arrows in Fig. 2d, and upper arrows 
in Fig. 4 and Fig. 5a. Since this feature has been already 
discussed with regard to the interaction of electron to 
sharp bosonic mode(s)^i^i^i^i^ii^i^iiS, we will not go into 
the details of this feature, except commenting upon its 
interesting doping dependence. For LEK, the size of this 
feature, which is interpreted as strength of electron-boson 
couphng, reduces upon doping in LSCO^ and Bi2201— 
while the size of HEA defined previously increases upon 
doping. It is then intriguing to ask whether an interplay 
of these two scales of low and high-energy anomalies will 
affect our understanding of the doping dependent effects 
seen in cuprates. 

IV. DISCUSSION 

The presence of three energy scales in the same data 
set hints at the hierarchy of interactions that are im- 
portant to the dynamics of electrons in cuprates. Aside 
from LEK which we believe is caused by electron-phonon 
interaction^iSiiSiiiS, HEA and the expanded band width 
are new observations that require more discussion. 

To gain more insights into the nature of the energy 
scales observed, Fig. 5 compares the data of the super- 
conducting (SC) sample, OP Bi2201, with that of the an- 
tiferromagnetic (AF) parent compound CCOC and ferro- 
magnetic (FM) Lai.2Sri.8Mn207 (LSMO). Comparison 
of LSMO and cuprates may give us some insight since 
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FIG. 5: Comparison of electronic structures of three per- 
ovskites in (a) superconducting (SC), (b) antiferromagnetic 
(AF) and (c) ferromagnetic (FM) phases. The images show 
ARPES data of (a) optimally-doped Bi2201, (b) undoped 
CCOC, and (c) LSMO. 

LSMO is reported to have similar pseudogap behavioral 
The LSMO comparison yields two insights: i) HE A is 
probably related to antiferromagnetism as it is absent 
in the ferromagnetic state and ii) the LDA calculation 
apears to give a correct band width if antiferromagnetism 
is not present. As for LEK near 0.03-0.09 eV, it is seen in 
metallic cuprates and LSMO but not in insulating CCOC 
where the polaron effect is too strong, thus suppressing 
the quasiparticle weight dramatically^^. 

HEA may be related to the short range Coulomb in- 
teraction as in the Hubbard (or t-J) model. Calculations 
using these models show that in undoped Sr2Cu02Cl2, 
the quasiparticle band width of the part below HEA is 
set by the J scale to be around 2-3 J ~ 0.25-0.35 eV^ 
while the higher energy part is presumed to be the inco- 
herent band with t scalc''^*. HEA in the Hubbard calcu- 
lation, which may come from the meeting of the quasi- 
particle band and the incoherent lower-Hubbard band, 
can be seen from the slightly underdoped to overdoped 
regimeSi; however, a further check still remains in the 
small doping regime where the calculation is challenging. 
In SrVOs, a feature similar to HEA is seen in ARPES^'"' 
and in LDA-I-DMFT calculation?'^,. With these models, 
the size of HEA is expected to be less pronounced upon 
dopingSl which is opposite to the doping effect we see in 
Fig. 4. 

One possible way to understand the anomalous dop- 
ing dependence could be the interplay between electron- 
electron and electron-phonon interactions. This drives 
the system into the polaronic regime in the underdoped 
samples2^, making a quantitative analysis of the size of 
the HEA difficult because obtaining polaronic physics 
correctly via an MDC analysis is challenging. Since the 
polaron physics is very strong in lightly-doped regime, 
this may artificially suppress the HEA size. 

Another possibility is that HEA comes from the renor- 
malization effect by bosonic mode(s). The in-plane plas- 
mon and the two-magnon mode are two possible can- 



didates. The coupling strength of in-plane plasmon in- 
creases upon doping, as shown by the energy-loss func- 
tion in Refi^. Qualitatively, the in-plane plasmon should 
give a very similar doping effect to that shown in Fig. 
4. The problem is that the plasmon mode energy in 
LSCO is around 0.8 eV, or twice as large. On the other 
hand, the two-magnon mode energy is right in the win- 
dow for the hole-doped cuprates ( 0.3-0.4eV) but not 
for F0234, electron-doped band ( 0.5-0.6 eV) while its 
strength quickly reduces upon dopin^S. 

Given the above alternative interpretations and issues 
related to them, despite the quantitative problem of the 
HEA size as a function of doping, the anomalous energy 
scale is still likely caused by the Mott-Hubbard physics. 
The fact that HEA is absent in FM LSMO (Fig. 5) 
but present in the Mott-Hubbard system SrVOa sug- 
gests that it is related to antiferromagnetism and Mott- 
Hubbard physics. 

Next is the enhancement of LDA band width which is 
anomalous because interactions usually reduce the band 
width. One possible candidate here could be the poor 
screening effect. With the Hartree-Fock equation of free 
electron, the Coulomb exchange term without the screen- 
ing effect^S will give rise to a negative term which results 
in an increased band width; the effect of poor screening 
can be seen in semiconductors. In cuprates, the electrons 
in lightly-doped systems are poorly screened compared to 
the overdoped regime and hence the discrepancy of band 
width is larger in underdoped systems. We note here 
that, for the SrVOs system, LDA-I-DMFT calculation^^ 
gives the band width correctly. This distinction between 
SrVOa and cuprates should be investigated further. 



V. CONCLUSION 

In conclusion, based on data from four families of 
cuprates over a wide doping range, we present evidence of 
a hierarchy of multiple energy scales in cuprates focused 
on: the low energy anomaly of 0.03-0.09 eV, a high energy 
anomaly of 0.3-0.5 eV, and an anomalous enhancement 
of LDA band width, extending over an energy scale of « 
2 eV. These results suggest that electron-phonon interac- 
tion, short-range Coulomb interaction and poor screen- 
ing should be all considered to understand the nature of 
cuprates. 
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APPENDIX: A FULL 2D ANALYSIS 

In the following, we will show that interestingly, a 
2D analysis could describe ARPES data well, covering 
a large energy and momentum space, by a compact set 
of parameters. This is an attempt to go beyond the con- 
ventional EDC or MDC analysis alone. However, we note 
that the physical meaning of such a parameterization still 
remains to be explored. 

Here we will apply the 2D analysis on the ARPES data 
of Pb-substituted Bi2201. The overdoped (OD) samples, 
Pbo.38Bii.74Sri.88Cu06+5, are non-superconducting (Tc 
< 4 K). ARPES data were collected with a photon energy 
of 42 eV. The energy resolution was set to 18 meV. The 
linear polarization of the light source is fixed to be in- 
plane along (0,0) to (7r,7r) throughout the measurement. 
Note that the fitted matrix element, which is shown in 
Fig. 6d-i, refers to this particular experimental geometry. 

The intensity measured in an ARPES experiment on 
a 2D material here will be parameterized by^^ 

/(k, LJ) = /o(k, i., A)f{uj)A{k, LO) (A.l) 

where Iq (k, A) is proportional to the one-electron ma- 
trix element and dependent on the polarization, momen- 
tum and energy of the incoming photon, f{Lo) is the Fermi 
function, ^(k, uj) is the single-particle spectral function 



given by 

[u-el- ReE(k, lo)^ + [ImE(k, ^ ' 

where is the bare band dispersion, and S](k, lj) is the 
self-energy. 

We note that we neglect the instrumental resolution 
here since the feature of interest is large compared to 
the resolution. In the following, we will assume weak 
momentum dependence of this extracted self-energy (i.e. 
S(k,c.) ^ E(c.) .) 

Fig. 6 shows the comparison of the raw ARPES data 
(first row. Fig. 6a-c) and the parameterized data (sec- 
ond row. Fig. 6d-e). We parameterize the ARPES data 
with the form given by Eq. Al where the spectral func- 
tion A{k^Ljj) is given by Eq. A2. In the fitting, we do 
not assume any form of the self-energy and the matrix 
element (i.e. every point of the fitted self-energy or ma- 
trix element is a free parameter in the fitting procedure.) 
The bare dispersion used here is the simple form given by 
the tight-binding (TB) parameters which is fitted to the 
LDA calculation shown in Fig. If and 2. The bare dis- 
persion is given by Eik) — —2t[cos(kxa) + cos(kyb)] — 
2t'cos{kxa)cos{kyb) — 2t"[cos{2kxa) + cos{2kyb)] — Ep 
where the TB parameters are t — 0.435, i' = — 0.1,<" = 
0.038, and Ep = —0.5231 eV. For simplicity, we use con- 
stant background in the fitting procedure. 

The fitting of the ARPES data (Fig. 6d-e) is very 
well in agreement with average error < 4%. Since the 
extracted self-energy does not show a strong momentum 
dependence, the approximate 3-dimcnsional spectra may 
be generated from this information as shown in Fig. 6i. 
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FIG. 6: On the first row, a)-c) are the experimental ARPES data along momentum direction as indicated by the band shown 
in bottom right. On the second row, d)-f) are the corresponding simulated images from the 2D analysis. The thicker color 
dashed lines are the dispersions generated from the tight-binding (TB) parameters given above. The solid color lines are the 
matrix elements (ME) obtained from the 2D analysis where the smaller white dashed lines are the empirically guessed form of 
the matrix element in the form of a cosine function (COS). On the third row, g)-i) are the MDCs of raw data (black dots) and 
corresponding simulated images (blue line), j) and k) are the extracted real and imaginary parts of the so-called self-energy in 
Eq. A2 for the cuts a, b and c; the extracted values are plotted in colors up to the energy not far from the bottom of the bare 
band (up to 0.6 eV for cut b and 0.25 eV for cut c) and in grey at higher energy. 1) shows the band structure generated from 
the TB parameters and the loosely-called self-energy along the nodal direction. 
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